UCLA/08/TEP/25 

Gravitational waves from fragmentation of a primordial scalar condensate into Q-balls 
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A generic consequence of supersymmetry is formation of a scalar condensate along the flat di- 
rections of the potential at the end of cosmological inflation. This condensate is usually unstable, 
and it can fragment into non-topological solitons, Q-balls. The gravitational waves produced by 
the fragmentation can be detected by Laser Interferometer Space Antenna (LISA), Advanced Laser 
Interferometer Gravitational- Wave Observatory (LIGO), and Big Bang Observer (BBO), which can 
offer an important window on the early universe and the physics at some very high energy scales. 



Supersymmetry is widely regarded as a likely candi- 
date for physics beyond the Standard Model. While 
many variants of supersymmetry have been considered, 
all of them have scalar potentials with some flat direc- 
tions lifted only by the supersymmetry-breaking terms. 
At the end of cosmological inflation, the formation of 
a scalar condensate along the flat directions can have 
a number of important consequences |l|. In particu- 
lar, it can be responsible for generation of the matter- 
antimatter asymmetry via Affleck-Dine (AD) mecha- 
nism [2J, and, in some models, dark matter can be pro- 
duced in the same process [3, [j]. Some flat directions 
could be responsible for the primordial inflation 5, 6]. 

The formation of AD condensate is a generic phe- 
nomenon, relying only on the assumptions of inflation 
and supersymmetry. In general, this condensate is un- 
stable: an initially homogeneous condensate can break 
up into lumps of the scalar field, called Q-balls [7j , under 
some very generic conditions [3(] . All phenomenologically 
acceptable supersymmetric generalizations of the Stan- 
dard Model admit Q-balls [||, which can be stable, or 
can decay into fermions p, Q • The formation of Q-balls 
is accompanied by a coherent motion of the scalar con- 
densate, which creates the source of gravity waves. We 
will show that fragmentation of the scalar condensate 
into Q-balls can produce gravitational waves detectable 
by LISA [13, LIGO III [Hj, and BBO [13]. 

The physics of AD condensate fra gme ntation has been 
studied both analytically bft fl3l [Tj, [lja, [lj| and numeri- 



cally [H H GJ, W, HI, Hi [231 ■ At the end of inflation 
(assuming that the inflation occurs in a hidden sector at 
high scales) the condensate has a uniform density, with 
small perturbations of the order of 10 -5 [l[. Under some 
rather generic conditions, the instabilities develop and 
lead, eventually, to the formation of Q-balls, which can 
either decay or remain as stable relics [3, [l3| . Although 
the final product of such evolution, Q-balls in the ground 
state, are spherically symmetric, the coherent motions 
associated with the condensate fragmentation and re- 
arrangement are not spherically symmetric. Moreover, 
the newly formed Q-balls first appear in their excited 
states and oscillate until theysettle in the spherically 
symmetric ground states [is Ell • The lack of spherical 
symmetry in the process of fragmentation is essential for 



generating the gravity waves. 

Following the general picture developed in Refs. [3lfla|. 
the scalar condensate undergoing fragmentation can be 
approximated, in the linear regime, as <p(x,t) = </>(t) = 
i?(£)e in C), plus a perturbation SR,SQ ex e s (*)~ lfcS . One 
finds that the homogeneous solution is unstable due to 
some exponentially growing modes, Re a > 0, where a = 
dS/dt 0, [la [lg|. The mass density of the condensate 
undergoing fragmentation can be written as 



p{x,t) =p +pi(x,t), 



where 



pi (x,t) = ep d k e ak cos(ujt - k ■ £ ) 



(1) 



(2) 



The instability develops when there is a band of growing 
modes with positive and large enough a,h [3J. The linear 
approximation breaks down when eexp(afei) ~ 1, but we 
will use this representation, up to its limit of applicability, 
to get the estimates of the gravity waves produced. 

The quadrupole moment that generates gravity waves 
is given by |24j 



D i:j = / d 3 x x l x 1 T m \ 



x,t). 



(3) 



where the energy-momentum tensor T 00 (x,t) fa p(x,t). 
The space integration is over some arbitrary volume. 

The power emitted in gravity waves in one frequency 
mode is given by: 



P(u) = \gu* fe»A» - ^|A»| 2 



(4) 



and the total energy emitted in gravitational waves, in 
all frequencies, is given by: 



E 



2n / p(uj)du> I x At 



(5) 



where At is the duration of the fragmentation. 

Based on the analytical and numerical calculations of 
the condensate fragmentation 0, [l8|, [25|, we take the 
typical parameters of the fastest-growing mode: 



£ fc x 10-.H*, oj k - vk ~ £ fe x 10 2 vH* 



(6) 



where H* is the Hubble constant at the time of the con- 
densate fragmentation, and v is the typical group velocity 
of the wave front in the evolution of the condensate, and 
we expect that the dimcnsionlcss factor ^ ~ 1, based on 
the results of Refs. [3j, Il8ll25j. 

Since no cancellations are expected in the absence 
of spherical symmetry, we replace the XiXj by (fk X 
10 2 il*) _2 in the space integration, take the volume to 
be V ~ H~ 3 , and assume that eexp(afet) ~ 1. Then, for 
the leading mode, 

Da (t) ~ H~ a (10 2 H„) ~ 2 p cos( Wfc i - kx) , (7) 



and, in frequency space, 



0„(uO~io-V|i 



(8) 



For uj ~ 10 2 vH„, we estimate the power in gravitational 
waves in a Hubble volume: 



P ~ 10 4 £ 



4,-2 Q 



m 



(9) 



To estimate the velocity of the wavefront in the process 
of fragmentation, we note that, for the mode (p(x,t) ~ 
R(t) exp{a fe t} cos(lj k t—kx), where R(t) is a slowly chang- 
ing function of time, 1 



■Sv 



£v a k /k, £ v ~ 1, 



(10) 



where the uncertainty factor £„ ~ 1 will be retained to 
keep track of the uncertainty in the final answer. 

The relation between a k and k is given by a dispersion 
relation [3(, which takes a simple form, 



(al + k 2 )(al-(n 2 -V"(R))) + 



m 



2„2 



= 0, (11) 



under the following assumptions 

n 



H < k ~ «fc < 
fi 2 — V"(i2) ~ tl ~ m^ valid in the case of the fastest- 
growing mode in gravity mediated supersymmetry break- 
ing models (the latter is essential for the (l ~ m^ con- 
dition 1]). This equation has an approximate solution: 
a fe k, fc/V3, and 



w fl ~£W*0 fl ~£(i/V3) ( 



io- 2 ^ ■ 



(12) 



The fragmentation takes place on the time scale of the 



order of At ~ a I 



C'io- 



'-ff* . (Here we neglect the 



possible contributions from the collisions and oscillations 
of Q-balls, which can take place on a much longer time 



1 The adiabatic limit R/R — > is amenable to perturbation the- 
ory Q. This case corresponds to a large global charge density. 
For a rapidly varying R(t), the numerical calculations give similar 
results regarding the fragmentation time and length scales [18H . 



scale 3,26].) The total energy in gravity waves generated 
in the Hubble volume is 



E-PAt-G^tfit 



Po_ £-3,6 

m 



(13) 



This corresponds to the energy density in gravitational 
waves at the time of production 



m-3f-3 C 6 Po 
Paw* ~ 1U £ fc $ v 



HIM 2 , ' 



(14) 



where Mp\ = 1/v8ttG is the reduced Planck mass. 
Hence, the fraction of the energy density in gravitational 
waves at the time of production is 



O 



GW* 



io- 3 C 3 ^ 



pI 



(H*Mpi) 4 



(15) 



If the energy density of the condensate is comparable 
to the total energy density, or if the condensate energy 
dominates the energy in the universe, then po ~ 3H 2 Mp V 
and flaw* ~ 10~ 3 . 

The energy density in the condensate depends on the 
model, and, foremost, on the type of supersymmetry 
breaking terms that lift the flat direction. This is be- 
cause the potential along the flat direction depends on 
supersymmetry breaking (it vanishes in the limit of exact 
supersymmetry), and there are many ways to break su- 
persymmetry. In gauge-mediated supersymmetry break- 
ing scenarios the potential can have the form [3j 



n.,),, u>„(i + M! 



(16) 



Here M$ is the scale of supersymmetry breaking, which is 
of the order of 0(1) TeV. In gravity mediated scenarios, 
the flat directions are lifted by mass terms that persist 
all the way to the Planck scale [13| : 

V{fp) « m% (l + i^log (j^r)) M 2 > ( 17 ) 

where K ~ 0.05 (for squark directions) describes the run- 
ning of the mass term [13| . Since m$ ~ Ms ~ 1 — 10 TeV 
in typical models, both potentials are phenomenologi- 
cally acceptable near the minimum. However, in AD 
condensate and inside the Q-balls that form in its frag- 
mentation, the vacuum expectation value (VEV) can be 
very large. 

The main difference between gauge and gravity medi- 
ated cases for us is the mass per baryon number stored in 
the AD condensate and in the Q-balls that form eventu- 
ally as a result of the fragmentation. In the gravity medi- 
ated scenarios, the mass density is po ~ rn^cj) 2 , the global 
charge density is uq ~ m$<fi, and the mass per unit 
global charge is of the order of m^, independent of the 
VEV 4>o. In gauge mediated scenarios, the mass density 
is po ~ wi, the global charge density is tiq ~ m$(f> 2 , and 
the mass per unit global charge is Po/uq ~ rrA/4> @jQ- 



The flat directions that carry a non-zero global charge 
Q = (B — L) contribute to the generation of the baryon 
asymmetry via AD process [2| ■ The requirement that 



r ]B 



ns/n^ 



10" 



implies that the total mass den- 



sity of such a condensate cannot be of the order of the 
total density of the universe in generic models. This is 
true in both gauge-mediated and gravity-mediated su- 
persymmetry breaking models. The gravity waves from 
the fragmentation of such a condensate are well below 
the capabilities of the current and planned detectors. 

On the other hand, there are flat directions whose 
baryon number B and lepton number L are equal to each 
other 0, 51 El- While the flat directions with B ^ L 
contribute to baryon asymmetry of the universe, those 
with B = L, have zero (B — L) density. Electroweak 
sphalerons destroy any primordial (B + L) asymmetry, 
and so the corresponding t/b+l = nB+L/n-y is not con- 
strained. It is possible that, at the time of the fragmen- 
tation, rjB+L 3* tjb- For B = L flat directions, there 
is no reason why p cannot be of the order of the total 
energy density. The fragmentation of such flat directions 
can produce a detectable level of gravitational waves. 

There are various examples in the literature of the flat 
directions that can dominate the energy density of the 
universe, while they do not contribute to (and are not 
constrained by) the observed baryon asymmetry of the 
universe. This is the case, for example, when the effec- 
tive mass for the phase direction is large during inflation, 
which results in the initial condition with a very small 
f2 for the scalar condensate 4>(x,t) = <j){t) = R(t)e lQ ^ . 
This is also the case when the inflation is driven by a 
flat direction, udd, LLe or NH U L [3, [fj, l3Q |- Another 
well-know example is a flat direction that acts as a cur- 
vaton and dominates the energy density of the universe 
at the time of oscillations and decay [3l[ . In all of these 
cases, the net global charge of the condensate is negli- 
gible, but the fragmentation can still occur and produce 
Q-balls and anti-Q-balls [25l. [27J] . 

Once the gravitational waves are created, they are de- 
coupled from the rest of the plasma. We can estimate the 
peak frequency of the gravitational radiation observed to- 
day, /* = uj k /2ir: 



I = 



J* — J* 
a 



a* 

a rh 



9s,o 



0.6 mHz 



£*& ( j'qq J 



5s,rh 

1/6 



1/3 






(18) 



r r i 



1 TeV 



10H, 



where we have assumed that a* « a r h, which also means 
that during the oscillations of the AD condensate the 
effect of the Hubble expansion is negligible. The val- 
ues of relativistic degrees of freedom are g s ,rh ~ 300, 
9s,a ~ 3.36. The subscript "rh" denotes the epoch of 
reheating and thermalization, while the subscript "0" 
refers to the present time. As we discussed, the typi- 
cal frequency of the oscillations of the AD condensate 
is uj k ~ 10 2 if» |, [H, i|. Then for T, h ~ 1 TeV 
(such a value of the reheat temperature is natural when 



the flat direction is responsible for reheating the uni- 
verse [30, |33), the frequency is of the order of mHz, 
which is in the right frequency range for LISA [lflj. A 
higher temperature T r h ~ 100 TeV corresponds to the 
LIGO III frequency range, 10 - 100 Hz [ll|. Signals in 
both of these ranges will be accessible to BBO [12]]. Since 
the supersymmetry breaking scale is related to the energy 
in the condensate, as well as the reheating temperature, 
LIGOIII and BBO could be in the position to probe su- 
persymmetry broken above 100 TeV, beyond the reach 
of Large Hadron Collider (LHC). 

The fraction of the critical energy density p c stored in 
the gravity waves today is 



O 



G\Y 



= Vow* — U^ 



» i \ / H, 

1.6V - 10- 5 /ioo x 1/3 



h 2 



n 



GW* 



(19) 



10- 8 £fc 3 £>- 2 



where clq and Hq are the present values of the scale fac- 
tor and the Hubble expansion rate. LISA can detect the 
gravitational waves down to flcwh 2 ~ 10~ n at mHz 
frequencies, while LIGO III is sensitive to ilcwh 2 ~ 
(10~ 5 -10" n ) in the (5-10 3 ) Hz frequency band. There- 
fore, the gravitational waves with flew* ~ 10~ 3 and 
QGwh 2 ~ 10 -8 from the fragmentation of an AD con- 
densate can be detected. The first results of our numeri- 
cal simulations (work in progress) appear to produce the 
gravitational wave signal that is somewhat weaker. We 
attribute the difference to the value of the uncertainty 
factor $, k 3 ^y, which is especially sensitive to the average 
wavefront velocity. 

As one can see from eq. (|14[) . the power generated in 
frequency w is proportional to 



c k 



-3 



10 2 #* 



Hence, the spectrum is strongly peaked near the longest 
wavelength, of the order of the Q-ball size 0, [l8|, |25| . 
The relatively narrow spectral width will help distinguish 
this signal from the gravity waves generated by infla- 
tion [33j, which are expected to have an approximately 
scale- invariant spectrum (and a smaller amplitude). Fu- 
ture numerical simulations will help refine the prediction 
for the signal from Q-ball formation, which can help dis- 
tinguish this source from a phase transition in the early 
universe [11] • LISA and LIGO III will be able to dis- 
criminate the gravity waves due to fragmentation from 
those of point sources, such as merging black holes and 
neutron stars, which have specific "chirp" properties [351 ] . 
Furthermore, the signal discussed here will not create a 
significant background for the cosmic microwave polar- 
ization experiments, such as B-Pol 36], which can detect 
the gravity waves with extremely long wavelength. 

Some additional gravitational waves can be generated 
by collisions and oscillations of Q-balls 0, [26|] . We leave 
the discussion of the magnitude of this additional contri- 
bution to future studies. 



To summarize, the fragmentation of a scalar conden- 
sate into Q-balls, which is a generic consequence of super- 
symmetry and inflation, can produce a detectable level of 
gravitational waves, up to flawh 2 ~ 10 -8 , near the peak 
frequency of BOO and cither LISA or LIGO, depending 
on the reheating temperature. Detection of the gravi- 
tational waves form this process can shed light on the 
earliest post-inflationary epoch in the history of the uni- 
verse, can probe supersymmetry even if it is broken at a 



scale above 100 TeV, and can provide information about 
new physics at some very high energy scales associated 
with the flat directions. 
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